Finite-size effect of antiferromagnetic transition and electronic structure in LiFeP04 
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The finite-size effect on the antiferromagnetic (AF) transition and electronic configuration of 
iron has been observed in LiFeP04. Determination of the scaling behavior of the AF transition 
temperature (Tjv) vs. the particle size dimension (L) in the critical regime, 1 — •^^f^^-f , 
reveals that the activation nature of the AF ordering strongly depends on the surface energy. In 
addition, the effective magnetic moment that reflects the electronic configuration of iron in LiFeP04 
is found sensitive to the particle size. An alternative structural view based on the polyatomic ion 
groups of P04^ is proposed. 
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FIG. 1. (color online) (a) Iron ion environment of LiFeP04 
with P04~ ligand formed weaker crystal field. The P04~ 
polyatomic ion is constructed by four cr-bond and one res- 
onating TT-bond as shown on the right, (b) Alternative view of 
LiFeP04 with lithium residing in the diffusion channels along 
b-direction of large pore size, which has van der Waals-like 
bonding with the neighboring Fe-P04 groups. 



LiFeP04 has been considered as an ideal cathode elec- 
trode material with a workable redox potential of 3.45 
V and a high theoretical capacity of 170 mAh/g.[l| 
LiFcP04 is a band insulator in bulk crystal form in agree- 
ment with the calculated large band gap,[3 however, 



the fast electron/ion charge/discharge process has been 
shown to be possible only when the particle size is re- 
duced to below ^100 nm.Q While significant progress 
has been made to demonstrate the impact of nano- 
scaling on defect chemistry, electrochemical behavior, 
electron-ion conductivity, lithium ion diffusion, and the 
charge/discharge rate in general, 0^ the fundamen- 
tal question of "why nano-scaling matters?" remains. 
Generally, there have been two major approaches to 
interpret the property improvement for LiFcP04 us- 
ing nano-scaling induced electron/ion transport; the 
first is called the "domino-cascade" model, which de- 
scribes a fast anisotropic lithium insertion/extraction 
action coupled to the LiFeP04-FeP04 interface move- 
ment through an electron polaronic hopping mechanism 
in nano-crystallites,[5| and the second focuses on the sur- 
face effect with carbon coating and miscibility gap size 
reduction as a result of nano-scaling. 0, 0] 

Conventionally, the iron ion has been viewed as Fe^+ 
with oxygen ligands to form the octahedral crystal elec- 
tric field (CEF) of eg-t2g splitting. Q However, the mag- 
nitude of the effective magnetic moment (/ie//) obtained 
from the Curie- Weiss law fitting was inconsistent and 
spread in the range of /^e//~4.7-5.5 /is, which has been 
interpreted that Fe^"^ must be in the high-spin state (HS) 
as tf^e^ of S=2 with a spin-only value of fJ-eff ~ 4.9/iB, 
and the inconsistent moment above 4.9 fiB has often been 
described by incomplete orbital-quenching, i.e., the the- 
oretically calculated value from the Lande g-factor with 
J=L-|-S (S=2 L=2) without orbital quenching should be 
Me//=6.7^B.[3, Considering a typical iron ion in the 
crystal field composed of oxygen ligands with high elec- 
tronegativity, the HS choice of d^ implies an unexpect- 
edly small crystal field splitting between Cg — t2g, un- 
less the influence of the nearby P04~ phosphate groups 
arc included. In addition, within the explanation of in- 
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complete orbital-quenching, the actual mechanism that 
triggers the different degree of orbital-quenching remains 
unknown. Several alternative but inconclusive explana- 
tions for the excess moment have been suggested, includ- 
ing the existence of the FeP04 phase, hidden impurity 
phases with Fe"^"*" , small magnetic polarons of large polar- 
ized macro-spin. |1 2] and the surface disorder-generated 
low-spin state (S=l/2) of Fe^+.fl^ 

In the discussion of the spin state of the iron ion, we 
find the traditional view starting from the distorted FeOg 
octahedral crystal electric field cannot avoid the influence 
of the nearby phosphate group of PO^". Alternatively, 
weaker crystal field surround Fe can be formed by the 
weak ligands of (P04)^~ as polyatomic ions as shown in 
Fig. [TJa). Owing to the extremely high positive charge 
density of hypcrvalent P^+ with a particularly small ion 
radius (0.17A (IV)), the (P04)'^~ group has often been 
treated as one polyatomic ion of tetrahedral coordina- 
tion in chemical reactions. The standard view of the 
phosphate group P04~ is shown in Fig. [TJa), which is 
composed of P-0 tetrahedral-shaped sp^d hybrid orbitals 
with four (T-bonds plus one resonating 7r-bond. 

Considering (P04)'^~ as a unit, the larger effective ra- 
dius leads to a much lower average negative charge den- 
sity per phosphate than pure 0^~. The actual environ- 
ment surrounding the iron ion can be described as either 
an edge-sharing Fe-(P04) pair, or an Fe-(P04)5 crystal 
field as shown in Fig. [1] While the neighboring oxygen 
atoms of high negative charge density for FeOg octahe- 
dron lead to strong CEF splitting of eg-t2g ^ 1-2 cV, the 
CEF constructed by the five phosphate groups of much 
lower electronic charge density and longer average Fe- 
(PO4) distance would create a weak crystal field, which 
concurs with the experimental finding of the HS state 
choice for Fe-^+ of d^. The ID diffusion channels along 
the b-direction of large pore size for lithium could pre- 
fer van der Waals-like bonding as depicted in Fig[ljb), 
similar to that of the typical battery cathode Lia;Co02 
with lithium in the van der Waals gap between 2D C0O2 
layers. 

A change at the electronic configuration level allows 
LiFeP04 to transform from a bulk band insulator to a 
good conductor through nano-scaling alone. To explore 
this mechanism using the newly proposed structural view 
based on the polyatomic ion unit of P0|~, we report a 
complete magnetic property study of LiFcP04 with con- 
trolled particle size, starting from a bulk single crystal to 
nano-scaled particles, of ~50 nm level, in discrete steps. 
The effective magnetic moment of iron has been found to 
depend on the particle size. The effective magnetic mo- 
ment changes as a result of the increase in surface energy, 
which implies that a subtle electronic structure change 
for iron occurs. This result could be closely related to 
the observed significant improvement of electronic con- 
duction for the cathode material LiFeP04, which is com- 
posed of nano-scaled grains. Q In addition, we find that 
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FIG. 2. (color online) (a) Average spin susceptibilities as a 
function of temperature for various particle sizes. The single 
crystal data were obtained from the powder average of the 
oriented measurements as described in the text. The isother- 
mal magnetizations M(H) at 300 K for the single crystal along 
three major orientations are shown in the inset, no sponta- 
neous moment from FM impurity is found near the zero field. 



the antiferromagnetic (AF) phase transition temperature 
T/v in the critical regime shows a scaling behavior in 
agreement with that predicted by the finite-size effect. 
The particle size dependence of the electronic activation 
type has been established rigorously for both the para- 
magnetic spin state and the critical phenomenon for AF 
ordering. 

A series of samples with different particle sizes (from 
bulk single crystal to ^50 nm nano-crystallites) was 
prepared from one crystal batch using the jet-milling 
method. The experimental details of crystal growth, 
particle size reduction method, effective magnetic mo- 
ment analysis, and particle size distribution are described 
elsewhere. (l3| Chemical analysis using the ICP method 
has confirmed the Li:Fe ratio to be 1.002(5):1. The aver- 
age spin susceptibilities (M/H) as a function of tempera- 
ture for various particle sizes are shown in Fig. [21 a). The 
effective magnetic moments (fJ-eff) were extracted from 
the Curie- Weiss law fitting in the paramagnetic regime 
(100-300K) as X = Xo+C/(T-e) directly with the Curie 
constant C = No/ig^^/Sks. Although Curie constant can 
also be extracted from the conventional approximation 
using l/x{T) vs. T linear fitting, the nonlinearity intro- 
duced by the no-longer-negligible Xo contribution at high 
temperature can be avoided. In particular, to maintain a 
strict particle size dependence analysis of fJ-eff, the bulk 
single crystal sample has not been crushed into a powder, 
and its Curie constant has been extracted directly from 
the anisotropic spin susceptibilities along three major or- 
thogonal axes following the equations derived by Liang 



3 



P X Specific surface area (cm^/cm^) 



10000 
A 




2.4 



2.3 



"^■^ CD 



- 2.1 



- 2.0 



0.1 1 10 

Particle size (\xm) 



100 



FIG. 3. (color online) Particle size dependence of the effec- 
tive magnetic moment fi^f / , where the /^e// values are derived 
from the fitted Curie constants using the Curie- Weiss law of 
M/H(T) shown in Fig. [S] The upper x-axis shows the corre- 
sponding specific surface area with an undetermined geomet- 
ric factor of P {P=6 for cubic). The corresponding g-factors 
are estimated by assuming Heff=4:.9fj,B of S=2. 
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The effective magnetic moment vs. particle size is sum- 
marized in Fig. [31 virhich has been analyzed from the av- 
erage spin susceptibility M/H(T) (Fig. (2) with an applied 
field of 0.5 Tesla. It is clearly demonstrated that ^e/ / for 
the bulk single crystal is ^^5.9 and this value mono- 
tonically reduces to ^--^4.9 for the nano-sizc (~50 nm) 
particles. The HS-like configuration for Fe^+ of d*' is 
configured as e^t^g with S=2 of ^e//=4.9 /is (assuming 
g=2). Because /Lte//=4.9 ^b and ^eff=(!i-7 f^B repre- 
sent the theoretically calculated values using J=S and 
J=L-|-S (with L=S=2), respectively, the figff values 
larger than 4.90 /is could tentatively be interpreted as 
a result of the incomplete orbital-quenching. Alterna- 
tively, we may infer the change in Heff from changing 
the g-factor beyond the electron g-factor of 2 as shown 
in Fig. [3l which could be described as a result of change 
in local field. A monotonic particle size dependence of 
/ie// has been constructed as shown above, which im- 
plies that the earlier reports of inconsistent /ie// values 
arc closely related to the particle size. However, the fol- 
lowing questions arise: why docs the perfect spin-only 
value of /ie//=4.9 jUs exist only in the nano-scaled par- 
ticles, and why is fJ-eff higher in the single crystal sam- 
ple? Above all, whether this phenomenon is described by 
incomplete orbital-quenching or local field change, there 
must be an intrinsic electronic structure change as a func- 
tion of particle size. 

As shown in the inset of Fig. [21 the isotherm M(H) data 
at 300K have completely ruled out the existence of fer- 
romagnetic impurity due to the absence of spontaneous 
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FIG. 4. (color online) The scaling behavior of reduced tem- 
perature vs. dimension is plotted in the log-log scale, and the 



critical exponent u for 1 
to 1 



Tjv(£) 



L is found to be close 
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For the large and nano-size particles, the high surface 
energy state is shown in red, and the ground state is shown in 
yellow. The inset shows the AF spin gap and its relationship 
to the surface energy, as described in the text. 

magnetization jump near the zero field. We find that the 
temperature-independent part of the average spin sus- 
ceptibilities (xo) is still at least one order higher than 
that estimated from the core diamagnetic, the Van Vleck 
paramagnetic, and the Pauli paramagnetic contributions 
combined, [l^ In addition, Xo level decreases monotoni- 
cally with the particle size as indicated in Fig. [2l which 
clearly opposes the assumption of higher itinerant car- 
rier densities for smaller particle size, i.e., a higher t^^°"'* 
contribution due to the increased carrier density of states 
near the Fermi level. These results strongly suggest that 
a possible spin state change may have happened following 
the particle size change. 

In addition to the particle size dependence of the effec- 
tive moment in the paramagnetic state, Tjv has shown an 
insignificant but consistent reduction for smaller particle 
sizes as revealed in Fig.[2{a), with T^ ^52 K for the sin- 
gle crystal and ^49 K for the nano particle of ^50 nm. 
We find that the relationship between Tn and the parti- 
cle size can be scaled nicely by the reduced temperature 
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particle size dimension L, where T;v(oo) 
corresponds to the single crystal sample of assumed infi- 
nite dimension, as 
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and ly is fitted to be near 1 as plotted in Fig. |4l Such 
scaling behavior perfectly demonstrates the finite-size 
effect when the magnetic correlation length is limited 
by the system dimension due to the effect of geometric 
confinement. [l7 . 18 1 The current extracted critical expo- 
nent v^l agrees well with the mean field theory prcdic- 
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tion of a 3D Heisenberg system at the critical regime. 

We find that the observed scahng behavior of the crit- 
ical exponent i/^l affects the role of surface energy in 
a nano system. Let us consider a single crystal sample 
of volume V, which is reduced into nano size particles of 
spherical shape with radius r directly through jet milling 
particle size reduction. We also assume that the observed 
reduction of Tp^{r) at nano size is due to the logarith- 
mically enhanced surface energy of jA (A=total surface 
area) from the bulk (7A^0) of a well-defined Tjv, i.e., 
it requires less thermal energy ksTiy^r) to cross the AF 
spin gap [ksTiq) in the critical regime, as illustrated in 
the inset of Fig. U) Then we can derive the scaling be- 
havior directly in the mean field approximation as 

ksTN - ksTNir) - 7^ = 7[( 4^)(47rr2)] ~ r-\ (2) 

which is comparable to the dimensionless scaling behav- 
ior of the critical exponent i'=l as shown above. In the 
critical regime of AF transition, the correlation length 
of the order parameter, i.e., the stagger magnetization 
of the AF ordering, would diverge logarithmically with 
the reduced temperature. Previously, the critical phe- 
nomenon was revealed by the scaling behavior using the 
magnetic correlation length (^). Now, this phenomenon 
can be explained intuitively using the role of surface en- 
ergy in the nano system, i.e., the increased surface energy 
plays an equally important role as the thermal energy in 
a nano system, which should in turn adopt a different 
electronic structure at a raised ground state higher than 
that of a bulk single crystal, as reflected clearly by the 
changing fief / values as shown in Fig. [3l 

In summary, by carefully analyzing the effective mag- 
netic moment of samples from single crystal to nano-size 
of '^SO nm, we find that the finite-size effect on the tran- 
sition temperature and the effective moment of LiFeP04 
reflect the significantly enhanced surface energy. The 
observed size dependence of fieff should be viewed as a 
reflection of the change in the electronic structure as a 
result of competing thermodynamic factors to lower the 
total Gibbs free energy, i.e., the surface energy 'fAA, 
the entropy TAS, and the magnetic energy HAM in 
AG = AU-TAS+jAA-HAM. Each term in the total 
Gibbs free energy should be given nearly equal weight at 
the nano scale to minimize AG in an intricate and com- 



petitive manner, instead of the enthalpy change dom- 
inating the minimization, which only reflects the ideal 
electronic structure. 
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